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Abstract 

Wc study the propagation of electromagnetic waves in a chiral fluid, where the molecules 
are described by a simplified version of the Kuhn coupled oscillator model. The eigenmodes 
of Maxwell's equations are circularly polarized waves. The application of a static magnetic 
field further leads to a magnetochiral term in the index of refraction of the fiuid, which is 
independent of the wave polarization. A similar result holds when absorption is taken into 
account. Interference experiments and photochemical reactions have recently demonstrated 
the existence of the magnetochiral term. The comparison with Faraday rotation in an achiral 
fiuid emphasizes the different symmetry properties of the two effects. 



Resume 

On etudie la propagation d'ondes clectromagnctiques dans un fluide chiral, dont les 
molecules sont decrites a I'aide d'une version simplifiee du modele d'oscillateurs couples 
de Kuhn. Les modes propres des equations de Maxwell sont des ondes polarisees circu- 
lairement. L'application d'un champ magnetique statique entraine I'existence d'un terme 
magnetochiral dans I'indice de refraction du fluide, independamment de la polarisation de 
I'onde. Un rcsultat semblable s'applique au cas de I'absorption. Le terme magnetochiral a 
ete recemment mis en evidence dans des experiences d'interferences et dans des reactions 
photochimiques. La comparaison a la rotation Faraday dans un fiuide achiral souligne les 
differences de symetrie des deux cas. 
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I. INTRODUCTION 



In simple fluids, the propagation of an electromagnetic wave is usually treated as fol- 
lows. Sticking to a classical description, one considers that the bound electrons of the fluid 
molecules are displaced from their equilibrium positions by the Lorentz force of the wave- 
fleld. This (time dependent) displacement or induced electric dipole can be modelled as a 
current density j(r, t) which is related to the electromagnetic field E(r, t) and B(r, t) through 
Newton's force equation (Complications such as local fields effects will not be considered in 
this paper). Neglecting the magnetic contribution (the Bohr model of the hydrogen atom 
leads to an estimate I I ~ l^^l ~ - < 10~^), and using Maxwell's equations then leads 

^ elec E ^ 

to textbook expressions for the index of refraction of the fluid [1] 

For chiral molecules, one has to consider the displacement of the bound electrons on a 
well defined geometric structure. In other words, one has to take into account the variation of 
the electromagnetic field down to molecular distances. The induced electric dipole moment 
(or current density) depends now on the field and on its spatial derivatives. This may be 
the reason why the propagation of electromagnetic waves in a chiral fluid is not frequently 
treated at an undergraduate level [2]. 

On the other hand, students may encounter phenomenological models which show that 
Maxwell's equations in an isotropic chiral fluid have circularly polarized waves as cigcnmodes 
[3-5]. This in turn suggests a comparison with the Faraday effect [6], which deals with the 
propagation of waves in an simple (achiral) fluid in the presence of a magnetic fleld. The 
comparison of these two situations raises questions about parity and time reversal transfor- 
mations in a "non particle physics" context [7-10]. These symmetry considerations are not 
confined to theory: recent experiments [11-13] have shown that the interplay of chirality and 
magnetism have remarkable consequences, that have been observed in interference experi- 
ments and in photochemical reactions. Roughly speaking, the application of a magnetic field 
on a chiral fluid leads to a change in the index of refraction, and this new "magnetochiral" 
term is independent of the state of polarization of the wave. The order of magnitude of this 
term is rather small (see below), but its possible implications for the origin of terrestrial 
homochirality in biomolecules (DNA, proteins,...) are rather interesting [14]. 

The present paper does not claim to be original. Its merit is perhaps to bring to the 
attention of students some topics which are somehow scattered through the litterature. 
Furthermore, experiments which are both "hot" and accessible to undergraduates are not 
so frequent. 

To illustrate these considerations in a self contained way, section II will present polariza- 
tion conventions and notations, as well as some deflnitions pertaining to chirality. Section 
III presents a sketchy derivation of the Faraday effect in a simple (achiral) fluid. In the 
framework of the clastically bound electron model, this model yields reasonable orders of 
magnitude, and enables one to study various aspects of the Faraday effet (dependence on 
the direction of propagation of light with respect to the magnetic field, refiection on a con- 
ducting mirror). The introduction of a damping term in the electron's equation of motion 
can then model absorption phenomena. 

In section IV, we consider a model of chiral molecules, which consists of two coupled 
anisotropic oscillators [15,16]. By further restricting the accessible orientations (albeit in an 
isotropic way), we will derive to lowest order the infiuence of chirality on Maxwell's equations 



3 



and show explicitly that the eigenmodes are circularly polarized waves. The differences with 
the Faraday effect will be pointed OTit, in particular from the symmetry point of view. 

In section V, we consider the effect of a magnetic field on a chiral fluid, along the lines 
of Larmor's theorem. This non-rigorous point of view [10] suggests the existence of a new 
(symmetry allowed) term in the index of refraction, called the magnetochiral term. Its 
order of magnitude is obtained through the model of section IV, and we discuss a recent 
interference experiment where this term is involved. Finally, when absorption is taken 
into account, the magnetochiral contribution leads to remarkable results in photochemical 
reactions. 



II. DEFINITIONS AND NOTATIONS 



Two molecules are superposable if one can bring them into coincidence by using only 
translations and rotations. A molecule (L) is chiral if it is non-superposable onto its mirror 
image (D). The (L) and (D) forms of a molecule are called enantiomers of this molecule. 
A fluid made of equal amounts of both enantiomers is called racemic. A fluid made only of 
(L) (resp. (D) ) molecules will be called a (L)-fluid (resp. a (D)-fluid) in this paper. 

In agreement with many authors [16-19], we consider a wave to be right (resp. left) 
circularly polarized wave if an observer looking towards the source sees at a flxed point in 
space, the tip of the electric field turn clockwise (resp. anticlockwise) with time. 

Physical properties pertaining to a right (resp. left) circularly polarized wave will be 
denoted by the subscript + (resp. -). A subscript L (resp. D) will be further added if one 
deals with a (L)-fluid (resp. (D)-fluid). So 11^(00) denotes the index of refraction of a right 
circularly polarized wave of frequency cu in the (L)-fluid. 



III. THE FARADAY EFFECT 

A. The elastically bound electron model 

We consider a monochromatic plane wave, of frequency cu, propagating in a simple (achi- 
ral) fluid. The direction of propagation is the 2; > direction. We model the atoms as 
hydrogen atoms, where the bound electrons are submitted to an elastic force / = —rrieUiQf, 
where rrie is the electron mass, luq a typical electronic frequency (a;o ~ 10^^ Hz), and r is 
measured from the equilibrium position of the electron. Neglecting the magnetic force, we 
get 

nier = —meU)Qr — eE (1) 
If N is the particle density, the induced electric dipole moment per unit volume then reads 



P^-7Vef=+ .^f^ ^. (2) 



leading to a polarisation current density j = Plugging these results in Maxwell's equa- 
tions [20] , and denoting by c the speed of light in vacuo yields 
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Looking for a plane wave solution E = Re(Eoe-^^'^^~'^*''), and introducing the index of 
refraction n{uj) through k = ^7^^ , we get 

n\u;) = l + 4 



B. Effect of a magnetic field and circular polarizations 

— * 

We specifically consider the propagation of a right circularly polarized wave (E+ = 
Eo(cos(A;2; — out), sm{kz — cut), 0)) in the presence of a static magnetic field Bq — BqCz- The 
equation of motion 

rrier — —meCUQr — e{Ej^ -\- r X Bq) (5) 

shows that the Lorentz force due to Bq is antiparallel to the elastic force. Following the 
same steps as before gives 

nl{u, So) = 1 + . (6) 

me£ol<^0 — ^ — ^"^J 

— * 

For a left circularly polarized wave (E_ = Eo(cos(fc2; — out), —sm{kz — cut), 0)), we find 



n^iu;, Bq)^1+ (7) 



since in this case, the Lorentz force due to Bq is parallel to the elastic force. 



C. Faraday rotation 



Let us now consider a linearly polarized wave propagating in a confined fiuid (0 < 2; < /) , 
in the presence of the magnetic field Bq = BqCz- Since a linearly polarized wave can be 
decomposed into two circularly polarized waves (+) and (— ) of equal amplitude but with 
different phase velocities {v± — c/n±{uu, Bq)), elementary calculations show that the direction 
of vibration of the electric field rotates between z — Q and 2; = Z by an amount a, with 

a = ^ An{uj, Bq) (8) 

where ^n{uj, Bq) = n+{cu, Bq) — n_(cj, Bq). This rotation is counterclockwise for an observer 
looking towards the source since v+ < as shown by equations (6,7). (Note that this 
simple model applies only to diamagnetic materials; paramagnetic materials require a more 
sophisticated treatment) . 
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Equations (6,7) also show that, as long as one deals with frequencies ul — one 

has 

n±{uj,Bo) c^n{uj±ijjL) (9) 

up to 0(— )^ terms. Since a; ~ ~ 10^^ Hz and even for Bq = 10 T, we only have 
ujl ~ 10^^ Hz, we will consider that Larmor theorem, as expressed by equation (9), holds. 
This allows us to rewrite the magnetic birefringence as 

An a;, Bq) ~ ^ ' 10 

me duj 

Typical orders of magnitude for a liquid such as CS2 are An(a;, Bq) ~ 10~^ for Bq = 1 T. 
The observation of the Faraday birefringence in an interference experiment is relatively easy 
since the Rayleigh and Michelson interferometers may detect index variations down to 10~^ 
[21]. 

The above expressions pertain to a wave propagating parallel to the magnetic field Bq. 
For propagation in direction u (with m = |), equations (9,10) read 

— * 

n±(a;,5o)~n(a;±^^) (11) 

and 

A / „v e(Bn-u)dn(u) 

An a;. Bo) ~ ^ ° ^ ,^ ^ 12 

rrie duj 

as can be checked for a wave propagating in the z < direction {Bq ■ u < 0). 

The effect of placing a perfectly conducting mirror z — I va. the original experiment 
(where the incident wave has u — e^) can be analyzed from different points of view. The 
simplest one is probably described in [9]: for an observer receiving the reflected wave, the 
incident wave from 2; = to 2 = Hs equivalent, in the (symmetric w.r.t the mirror) image 
space, to a wave propagating in the z <{) direction in the same magnetic field Bq = BqCz (this 
is a clear illustration of the axial character of the magnetic field Bq, since Bq is perpendicular 
to the mirror). Taking into account the reflected wave, shows that Faraday's rotation is 
doubled for this observer (the reflection on the mirror is irrelevant since it changes the 
sense, but not the direction of the electric field). Typical orders of magnitude for a liquid 
such as CS2 are a ~ 0.2 - 0.3 rd, for / = 10"^ m and An(a;, Bq) ~ 10"^. 

Finally, the inclusion of damping in equations (1,5) leads to an absorption of the wave 
by the fluid. Denoting by 77,2 the imaginary part of the index of refraction, we clearly 
have n2+{uj, Bq) 7^ n2-{uj, Bq). This magnetic dichroism implies an elliptical polarization 
(together with a Faraday rotation of the major axis of the ellipse) in the above experiments. 
Typical orders of magnitude are '^"^^^^^"^ ~ 10^"^ — 10^^ for Bq = 1 T. 
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IV. PROPAGATION IN A CHIRAL FLUID 



A. The simplified Kuhn model 

As previously mentionned, to take chirality into account requires a rather detailed geo- 
metric description of the fluid molecules. We consider here the simplest model of a chiral 
molecule [15,16], which consists of two coupled oscillators. 

To fix notations, we consider a fixed trihedron (Oxyz), with unit vectors (e^;, e^,, e^), 
and a (L)-fiuid with N molecules per unit volume. Each molecule has two electrons, whose 
equilibrium positions are called and We denote hy d — \E^ — Rl\ the size of the 
molecule. We also define ^R^-Rl^d &?2 and ^0 = 

Due to their interaction with the propagating wave, the electrons are displaced from R^ 
(resp. R^) by an amount fi (resp. r2). 

— * 

For a fixed orientation of the molecule (i.e. for a fixed 6°2) ) the displacements of the 
Kuhn model are both one dimensional and coupled: 

(i) the unit vector in the direction of rl (resp. is denoted by bi (resp. 62)- The 
chirality stems from the fact that the unit vectors 61, 62 et 6^2 are such that 

X = 6?2 • (h X h) = -1 (13) 

(ii) the potential energy of the two electrons is written as 

Ep = ^mecoliff + rf ) + me^^l^i ' ri){h ■ ^2) (14) 

where ujq and fii2 are electronic frequencies. 

Compared to the original Kuhn model, wc further add two restrictions [16]: 

(iii) With respect to the fixed (Oxyz) trihedron, each molecule of the fluid can only 
adopt the six orientations (Ai, i = 1, 2, 6) shown in Figure 1. 

The geometrical parameters of (Ai) are denoted with an index (i): we thus have 
bi\ 62 \ r^\'f^\ The chirality is clearly the same for all orientations 

X« = 6?!^ • (6? X 6^) = X = -1- (15) 

(iv) We assume that the orientations (Ai) are equiprobable. 

These restrictions preserve the isotropy of the fluid, and make the calculations easier 
(more complete calculations can be found in ref [16]). 

B. Equations of motion 

Neglecting the magnetic contribution to the Lorentz force, we get 

rriefi = -rrieUjlfi - mefll2bi{b2 ■ - eE(ri) (16) 

and 
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mef2 = -meu;lf2 - me^l2^2(bi ■ fi) - eE(f2) (17) 
Projecting on bi and 62, and looking for forced solutions, we obtain for orientation (Ai) 



J2, 



where the orientation independent matrix M is given by 



M 



''12 

02 / ,2 , ,2 



and M"^ is its inverse. 

Let us consider for the time being a linearly polarized wave, with E = Re(E(^) e~^'^*)ex- 
The resolution of the equations for each (equiprobablc) orientation yields the induced electric 
dipole moment P of the fluid. It will turn out that P has components both parallel (P//) 
and perpendicular (P_l) to the fleld E. It is easily seen that orientations {A^,Aq) do not 
contribute to P, and that orientations (A3,A4) only contribute to P//. On the other hand, 
orientations (Ai, A2) contribute to both components. Elementary algebra show that 

— * 

Note that the perpendicular dipole component Pj^ does not vanish because of the finite- 
ness of d (see Figure 1). To lowest order in d (or more appropriately in |), one gets 

where the gradient term is calculated at Rq. 

Gathering all contributions to P and considering a general state of polarisation of the 
wave, we get 

P = a{uj) E + 7^(0;) curl E (21) 

— * 

where the electric field and its derivatives are calculated at Rq. In equation (21), one has 
and 

. Ne' QlM-ihxh)) .... 

where we have used equation (13). Equation (23) explicitly shows that ^^{oj) is a pseudo 
(or axial) scalar. This feature of '^^{uj) is indeed required by equation (21), since P is a true 
(or polar) vector and curl E a pseudo (or axial) vector [23] . 
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C. Rotatory power 



Using equation (21) together with j = one finds 

We are now back to the traditional study of Maxwell's equations in a chiral fluid [3-5] . 
For circularly polarized waves in the long wavelength approximation (| <^ 1), we get from 
equation (24) 

niiu;) ^ niu) ^ (25) 

with n'^(uj) = 1 + ^ and A;o = -. 

A rough order of magnitude for the natural birefringence {An^{uj) = n^{uj) — n^(u;)) of 
a chiral liquid can be obtained from d = 10 A, iV ~ 5 10^^ m^^, tu = = ^ = tt 10^^ Hz, 
leading to An^{uj) ~ 3 • 10"^ 

Several points can easily be checked on equations (24,25) 

(i) the index of refraction is the same for a (e.g. right) circularly polarized wave propa- 
gating in the (L)-fluid, in the 2; > and 2; < directions. 

(n) for the (D)-fluid, one has n^{u;) — n^{u;), since ^^{co) — —7^(0;). 

(iii) if a linearly polarized wave propagates in a (L)-fluid, confined between z — and 
z — I, the direction of vibration of the electric field will rotate by an angle (3 given by 

/3 = ^An^c.) = ^^^M^Ii (26) 
2 2£q 

In our model of a (L)-fluid, the rotation is counterclockwise, and a typical order of magnitude 

is /9 ~ 1.5 rd for / --^ 10"^ m. 

(iv) if one places a perfectly conducting plane at 2; = /, we again follow [9]: for an 
observer receiving the reflected wave, the incident wave is equivalent to a wave propagating 
in the (symmmetric w.r.t the mirror) image space in the 2; < direction. In the image 
space, the image of the (L)-fluid is the (D)-fluid (again a clear illustration of the pseudo- 
scalar character of chirality). The full experiment for this observer amounts to a propagation 
in the z < direction, flrst through the (D)-fluid (incident wave), and then through the 
(L)-fluid (reflected wave), leading to a cancellation of the angle of rotation. Once again, 
the reflection on the mirror is irrelevant since it changes the sense, but not the direction of 
the electric field. This cancellation contrasts with the Faraday result and emphasizes the 
symmetry differences between the two cases [2,8,9] 

Finally, absorption of the wave can be modelled through damping terms in equations 
(16,17). Detailed calculations can be found in [16]. We only stress here that the symmmetry 
properties of the real and imaginary parts of the index of refraction are very similar. In 
particular, we have n2_|_(a;) 7^ n2_{uj) (natural dichroism) and n2±(a;) = n^^^luj). 
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V. EFFECT OF A MAGNETIC FIELD: THE MAGNETOCHIRAL EFFECT 



In principle, it is possible to study analytically the effect of a static magnetic Bq on equa- 
tions (16,17). It is much quicker, following [10] to apply Larmor theorem (11) to equation 
(25) and to write 

nl{u;,Bo)r^nl{u;±^^) (27) 
for a wave propagating in direction u. Expanding the right hand side of (27) leads to 

/, . 5 N _ ,N I ■ u) dnju) 7^(^)fco e d-f^{u) 

n^[uj, Bq) — n[u) ± — — =F — 7^ -j 1 [kq ■ Bq) [Zii) 

2me aw 2eo 'imeSo aw 

where ko — k^u = ^u. Note that the Larmor frequency shift has not be done on the "/cq" 
term of equation (25). This seems plausible since the magnetic field Bq will act only on 
the coupled oscillators, that is on 7(0;); the "/cq" term on the other hand, comes from the 
operator "curl" in equation (21) which is not affected by the dynamics. We refer the reader 
to reference [10] for further discussions on this point. 

The remarkable prediction of equation (28) is that the last term, called the magnetochiral 
term, does not depend on the polarization of the wave. It exists for linearly polarized or 
non polarized light, and has several interesting consequences: 

(i) magnetochiral birefringence: a (L)-fluid in a magnetic field and the associate (D)-fiuid 
in the same magnetic field do not have the same index of refraction for waves propagating 
parallel to the field, regardless of the polarization of the waves (one may also say that a (L)- 
fluid in a magnetic field does not have the same index of refraction for waves propagating 
parallel or antiparallel to the field). 

(ii) magnetochiral dichroism: the introduction of damping in the Kuhn model shows that 
the absorption coefficient of a (L)-fiuid in a magnetic field is not the same as the absorption 
coefficient of a (D)-fiuid, regardless of the polarization of the wave (one may also say that a 
(L)-fluid in a magnetic field does not absorb in the same way waves propagating parallel or 
antiparallel to the field). 

(iii) these magnetochiral effects vanish if light propagates perpendicularly to the field. 

(iv) For the mirror experiment described in section IV C, (iv), in the presence of a 
magnetic field Bq = BqCz, the contributions of the magnetochiral term cancel. The argument 
is the same as for the chiral contributions. 

From the Kuhn model, one may get an order of magnitude for the magnetochiral term 
of equation (28): setting d = lOA, = 5 lO^^m.-^, cj = = ^ = tt 10^^ Hz and Bq = 10 
T, and using equation (23), we find a contribution to the index of refraction of order 10"*^, 
quite close to the experimental limits [21]. One may further argue that this is precisely 
the order of magnitude of the terms which have been neglected in Larmor "theorem" (see 
equation (9)). So, one has finally to turn to experiments, which indeed have demonstrated 
the existence of the magnetochiral term(s). 
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A. Magnetochiral birefringence: an interference experiment 

Imagine a Young double slit experiment with linearly polarized light, and let Oz be the 
optical axis of the set-up. After slit 5*1, one adds: 

• a vessel Ci containing the (L)-fluid in a magnetic field Bq = BqCz followed by 

• another vessel Ci containing the associate (D)-fluid in the opposite magnetic field 

One adds after sht 5*2 : 

• a vessel C2 containing the associate (D)-fluid in a field Bq, followed by 

• another vessel C2 containing the (L)-fluid in the opposite magnetic field (— -Bo)- 

The vessels are identical. Analyzing this experiment with equation (28) shows that 
the phase shift of the interference pattern caused by the enantiomers is entirely due to 
the magnetochiral term. All other contributions cancel, whereas the magnetochiral term is 
multiplied by four. The experiment has actually been done with a Michelson interferometer 
in a slightly modified way [12]. 



B. Magnetochiral dichroism: a photochemical reaction 

Discriminating enantiomers is a major problem in chemistry. A classical experiment is 
the following (see [24,25] and references therein). One considers an initially racemic solution, 
that one illuminates with right circularly polarized light. The absorption of light gives rise 
to chemical reactions (D L) and (L D) [26]. In a certain range of parameters, these 
reactions are first order and we denote by ki and k2 their respective rates. Experimentally, 
one finds ki — X/q^F+I'^) ^2 = -f^/o^2+('^) where X is a constant and Iq the intensity 
of the light. The concentrations of the enantiomers [L] and [D] are given by 

ffl = -ffl = MD] - fe[Ll (29) 

At equilibrium, the solution is not racemic anymore, since one has 

/ [L] - [D] \ ^ ki-k2 ^ n^+{uj) - n\^{uj) ^ n\_{uj) - n\^{u) 

UL] + [D]^^ k^ + k2 niV(c^)+4+M n\_{uj)+n\+{uj) ^ ' 

so that the least absorbing enantiomer is in excess. 

This experiment has recently been done with natural (unpolarized) light, in the presence 
of a magnetic field parallel to the direction of propagation of light [13]: the magnetochiral 
effect implies that the absorption is not the same for the enantiomers, leading again to a non 
racemic solution at equilibrium. The experiment of [13], achieves a value y ~ 10~^, with a 
magnetic field Bq — IOT. 
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VI. CONCLUSION 



We have presented in a simple way old theoretical models and recent experiments, on 
"chirality, hght and magnetism" [27-30]. These topics can be of interest for undergraduate 
students, either from a physical or a (bio)chemical point of view, the more so since they 
also show up in other fields [31-34]. At a more advanced level, local field effects as well as 
quantum calculations can be introduced [35,36]. 

It is a pleasure to thank J- J. Girerd for discussions. 
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Figure Caption 



Figure 1: The allowed orientations (Aj) for a chiral molecule in the restricted Kuhn 

— * — * 

model. The displacement of the electrons are along bi and 62- We consider a hnearly 
polarized wave propagating along the z axis, with E = Re(E(z) e~^'^)ex, see text. In 

— * 

orientations (Ai) and (A2), where = e^, the electric field is not the same on the two 
electrons, leading to a non zero contribution to P_l, see equation (19). In orientations (A3) 
and (A4), where 6?2 = ^y? the electric field is the same on the two electrons, leading to a 
contribution to P//. In orientations (A5) and (Aq), where 612 = ^x, the electrons are not 
coupled to the electric field since (E • 61 = E • 62 = 0). 
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